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Abstract 

The well-known problem of unidirectional plane flow of a fluid in a half-space due to the impulsive motion of the plate 
it rests upon is discussed in the context of the second-grade and the Oldroyd-B non-Newtonian fluids. The governing 
equations are derived from the conservation laws of mass and momentum and three correct known representations of 
their exact solutions given. Common mistakes made in the literature are identified. Simple numerical schemes that 
corroborate the analytical solutions are constructed. 
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1. Introduction 

Following Truesdell and Rajagopalld2000l) . consider a fluid of density g(x, t) and velocity field u(x, t), where x is 



the spatial coordinate and t the temporal one. Conservation of mass dictates that 

g + g div u — 0, ( 1 ) 

where a superimposed dot denotes the material derivative: ^ = d$/dt + (u ■ grad)^. Cauchy's first law of continuum 
mechanics supplies the additional conservation of momentum equation: 

gu=divT + b, (2) 

where T(x, t) is the stress tensor, and b(x, f) represents the body force(s). 

We suppose that the fluid is incompressible and homogeneous so that g = gradp = => g = const. =: go(> 0) 
for all x and t. Then, Eq. (Q]) implies that div u — 0, meaning that such a fluid can only undergo isochoric motions. 
Assuming a Cartesian coordinate system, where x — xi + yj + zk with unit vectors in the three coordinate directions 
i, j and k, one such motion is the unidirectional plane flow u = u(y, t)l, which clearly satisfies divM = 0. Finally, we 
assume no external forces act on the fluid: b = 0. All this means that the fluid fills the half-space y > with a solid 
plate lying in the x-z plane (i.e., at y = 0). The motion is uniform (translation-invariant) in the x and z directions. 

In 1851, Stokes considered a specific case of such a unidirectional plane flowQ He was interested in the case 
wherein the plate at y = is set into motion suddenly at time t — + . In other words, the plate's velocity is given 
by U(t) = U(t)H(t), where H(-) denotes the Heaviside unit step function, and U(t) is some smooth function, i.e., it 
possesses as many continuous derivatives with respect to t o n (— oo +oo) a s needed, that we are free to specify. Stokes 
himself made the distinction between U(t) and U(t) clear ( StokesL 1851 , p. 101), yet a bewildering array of papers 



from the 1990s and 2000s fail to take this into account. Thanks to the no-slip boundary condition, the fluid near the 
boundary assumes the velocity of the the plate, i.e., m(0, t) = U{t)H(t). Consequently, if one needs to compute the 
acceleration of the fluid at the plate (e.g., when applying the Fourier sine transform to a mixed derivative), the correct 
expression is 'jjj(0, t) = U'(t)H(t) + U(t)6(t), where <5(-) denotes the Dirac delta function. Of course, such equalities 
are meant in the sense of distrib utions as the generalized functions H{t) and 5{t) fail to have point-values everywhere 



(IKolmogorov and Fominll 19751 §21). An ubiquitous but inexcusable mistake is to drop the U(t)6(t) term. 
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1.1. Second-grade fluid 

The Rivlin-Ericksen fluids, also known as fluids of grade n dTruesdell and Raiagopal , I2000I Chap. 6), are a model 
of isotropic simple fluids of the differential type. Their constitutive relation can be written as an expansion in terms of 
the Rivlin-Ericksen tensors A*. For the incompressible second-grade fluid it takes the form 



-pi + S, 



S = yU Ai + a\Ao + o^A 



trD = 0, 



(3) 



where p is the isotropic (indeterminate) stress, I is the identity tensor, S is the extra (determinate) stress, A] = 2D, 
Afc+j = At + At gradM + (gradM) T Ajt (t>l),Ds| [gradii + (grad«) T ] is the symmetric part of the velocity gradient 
(sometimes referred to as the infinitesimal rate of strain) tensor, and a T superscript denotes the transpose. The 
constant yUo(> 0) is understood in the usual sense of fluid viscosity from Navier-Stokes theory, and the second-grade 
(co nstan t ) param eters a\ and a 2 will be discussed shortly. 

Tind ( 1963 ) was the first to consider plane flows of the second-grade fluid and give correct solutions using the 



Laplace transform and the Bromwich integral inversion formula. However, he did not consider Stokes' first problem. 
Following his derivation, for which he is rarely given credit, the stress tensor for our unidirectional plane flow has the 
component representation 

-p + a 2 {%f ^o| + aig 



[T] = 



- P + (a 1+ 2a 2 )(%) 



dtdy 







Substituting Eq. into Eq. ©, we obtain 



£><) 



du dp d 2 u 



dt 



d 3 u 

dx ^° dy 2 1 dydtdy ' 











dp „ du d 2 u 

oy oy oy 



= —f- 



dp 



(4) 



(5) 



Since the plate is infinite, translational invariance in the x-z plane implies that the pressure cannot depend on x or z 
(IFetter and Walecka, I2003L §61), i.e., dp/dx = dp/dz - 0. This leads to p — p(y, t) — p M + (a\ + 2ai){duldyY from 
the second equation in Eq. (O, where is the ambient pressure and at most a function of t. The determination of the 
pressure is always overlooked in papers on this topic, yet it is a fundamental part of the solution to this problem. 

Finally, we note the following thermodynamic restrictions: a\ > a nd a\ + a% — jPunn and Raiagopallll995b . 
When a\ < 0, the problem becomes ill-posed in the sense of Hadamard ( Coleman et al. , fl965l) . 



1.2. Oldroyd-B fluid 

Oldrovdld 19501) proposed a number of constitutive relations for incompressible fluids with fading strain memory 



(retardation) exhibiting stress relaxation. The so-called incompressible Oldroyd-B fluid is the one with 

T = -pI + S, S + iiS=^ (Ai+^ 2 Ai), trD = 0, 
where the upper-convected time derivative ( Oldrovd , 1950i Sec. 3(a)) is given by 

$ = $ - (grad u) d~d grad u + (div «)£. 



(6) 



(7) 



Here, Ai and are the relaxation time and retardation time, respectively, and po is the fluid's viscosity (as before). 
Noting that S = S(y, t) due to translation-invariance in the x-z plane, the second equation in Eq. © has the component 
form: 



J xx + A\ — A\d yx -T^ - A.\b xy ^ b xv + A-i M^yy-dy- o 



'dy 

Syx + Al-gf ~ MSyy£, 
Szx + M-gf -hSzy^ 



ilS,. 



Syy + Al 

S zy + A i-df 



dS r ~ 

I— 



dS v 



-w 2 (|) 2 

-po 




,, du , j d 2 u 








(8) 



2 



From the assumptions that prior to start-up the fluid is at rest, we have that S(y, 0) = 0, whence the equations for the 
components S xz , S zx , S vy , S yz , S zy and S a give 



Szx 



S yy = S yz = S zy = S zz = Vf > 0, y > 0. 



(9) 



Since the stress tensor must be symmetric, i.e., T = T T (=> S xy = S yx , in particular), the remaining equations are 



S xy + A l 



dSxy du d 2 u 
— =Mo - +fM) A 2 —, 



at 



d u r. idu\ 
- 2A x S xy — = -2p A 2 I — I 



(10) 



This derivation was first given by_ Tanner! ( 1962k . though many authors employ a clumsy and abbreviated version of it 
without giving him any credit. IOldrovdl(ll950T Sec. 4) presents a similar derivation for a problem in polar coordinates. 
Substituting the first equation in Eq. © into Eq. and recalling that S = S(y, t), we obtain 



Qo 



du 
di 



5p + dS xy 
dx dy 



Q _ dp | dSyy 
dy dy 



(11) 



As in Sec. II. II translation invariance in the x-z plane implies dp/dx = dp/dz = 0. Then, p = p(y, t) — p^ + S yy (y, t) = 
Poo thanks to Eq. (0 and the second equation in Eq. (fTTT i. Next, we can eliminate S xy from the first equation in 
Eqs. ( fTOb and ( fTTb by taking the y derivative of the former and the t derivative of the latter. Thus, we arrive at 



du 



d 2 u 



d 2 u 



dy 



dy 2 



d 5 u 
' dydtdy 



(12) 



The two non-trivial components S xx and S xy of the determinate stress can be calculated from Eq. (fTUl i once u is found 

from Eq. (fT2l . This completes the formulat ion of the problem. 

Considerations from thermodynamics ( Rajagopal and Srinivasal 2000l) restrict the values of the relaxation and 
retardation times to be such that < A\, though here we present some solutions also valid for Ao > A\. Causality 
requires that A \ > 0. Then, for the problem to be well -posed in the sense of Hadamard, a necessary (but not sufficient) 



condition is th at A 2 > 0. Experimental observations support all of these restrictions dToms and Strawbridget 11953 
01drovdlfl958h. 



2. Exact solutions by integral transform methods 



2.1. Second-grade fluid 

Defining v := po/go and a := ai/go and supplying Eq. (0 with the boundary condition discussed in Sec.Q]and a 
proper decay condition as y — > oo, we have the following initial-boundary-value problem (IB VP): 



du 
di 

u(0, t) = U H(t), 
u(y, 0) = 0, 



d 2 u d 3 u 
dy 2 dydtdy ' 



as y — > oo, 



(y, t) € (0, oo) x (0, oo); 

t > 0; 
y >0. 



(13a) 

(13b) 
(13c) 



Second-grade solution representation 1 (IChristov and ChristovL lin press!) . Assume a > 0. Using first the Fourier 
sine transform in y and solving the resulting ordinary differential equation in t with the Laplace transform, one obtains 



u(y,t) = U H(t) 



2 r 

k Jo 



sin(gy ) 



■ exp 



-vft 
l+af 



2a r°°^sm(6 
n J 1 + a, 



-v£t_ 

1 +af 



dt 



(14) 



While many authors have attempted to obtain this solution, Christov and Christovl ( in press ) show that all of them 
make the mistake of dropping the U(t)5(t) term (recall the discussion in Sec . [D in th e expr e ssion for the plate's 
acceleration when applying the Fourier sine transform. A series of papers dErdogarl 120031 : lErdogan and Imrak , 
20051 l2007allbh promulgates this error, while incorrectly claiming that there is a "deeper" mathematical reason that 
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their erroneous solution does not agree with the correct Laplace transform sol ution. Meanwhile, three recent pa 



pers ( IZierep and Fetecaul 120071 IZierep et al.l 120071 : IZierep and BohninpL 120081) make use of the incorrect version 



of Eq. j[14| to perform some manipulations renderin g their results erroneous. The transform error is also commit- 
ted in (Fe tecau and Fetecaul |2002; ShenetaL . l2006h . wher ein an incorrect so lution of Stokes' first problem for a 
second-grade fluid over a heated plate is obtained. Similarly. iKhan et al.l d2008l) obtain erroneous solutions for MHD 
second-grade fluid flows. 

Second-grade solution representation 2 ( PuriL 1984 ). Assume a > 0. Using the Laplace transform in t and the 
Bromwich integral inversion formula, one obtains 



u(y,t) = U H(t) 



n Jo \ V» V a -1 1 V 



(15) 



Note that Eq. (TBl i can be transformed into Eq. (T5[ by the sub stitution 77 = £ 2 /(l +cr$ 1 ). This solution representa- 
tion is correctly generalized to the case of a porous half-space bv lJordan and Puril d2003h. while an erroneous version 
of the porous-half-space problem for the related Burgers fluid is corrected in djordanll2010l) . 



Second-grade solution representation 3 (Bandelli, Raiagopal and GaldiL l 19951) . Assume a > 0. Using the Laplace 
transform in t, the quotient splitting technique of MorrisoA il956 ) and standard Laplace inversion tables, one obtains 



u(y,t) = U H(t) 



,-*«jr e^(2V^)erfc(^-)d, 



(16) 



where erfc(-) is the complementary error function and Iq(-) the modified Bessel function of the first kind of order zero. 

Note that both Bandelli et al. ( 1995 ) and Puril ( 1984 ) omit that the H(t) pre-factor. This is not a triviality because, 
as the following theorem shows, lim / _ > o+ u(y, t) + - u(y, 0) for this problem. This is kn own as the start-up jump and 
its physical significance is gracefully explained in ( Jordan and Puri , 2003 : Jordan , 2010l) . 



Theorem 1 dBandelli. Raiagopal and Galdilll995l) . Let u(y, t) be a function for which du/dy, du/dt and d 2 u/dtdy are 
all integrable at (y, t) = (0, 0), then \u(e, 0) — m(0, e)| — * as e — > + . 

The contrapositive of this theorem states that if \u(e, 0) - w(0, e)| ■/-> as e — > + (i.e., the initial and boundary data 
are incompatible), the solution to the IB VP may fail to have integrable either first and/or mixed second derivatives 
at (y, t) = (0, 0). This means that the solution itself is ill-behave d (singular) there. Ind eed, for Stokes' first problem 



we have incompatible initial and boundary data. Therefore, from (Bandell i et al.L 1 1 9951 Eq. (3.10)), we have \u(y, t) 



u(y, f)0 This fact about limits of singular functions is 



u(y, 0)| as t — > + , whence lim,^o+ u(y, t) + lim,^o- 
completely unrelated to integral transform methods. In the recent literature, one can find baffling statements such as 
"it was shown that the previous attempts to solve the problem by using the Laplac e transform tech nique are erroneous 
and that the method of the Laplace transform does not work for this problem" dErdoganl 120031) . However, such a 
statement is patently false as the two Laplace and one Fourier sine transform solutions above are all correct, they 
satisfy all conditions imposed, and they are equivalent. 

2.2. Oldroyd-B fluid 

Again, defining v := juo/f?o and supplying Eq. ( fT2l with the boundary condition discussed in Sec. Q] and a proper 
decay condition as y — > 00, we have the following IBVP: 



du du 2 



d 2 u ^ d 3 u 
dy 2 dydtdy ' 



u(Q, f) = U H(t), 
u(y, 0) = 0, 



it 



as y — > 00, 



(y, f) £ (0, 00) x (0, 00); (17a) 
f>0; (17b) 
y>0; (17c) 



"It appears that lTannerHl962l) was the first one to realize this for the Oldroyd-B fluid. He notes tha t "the integra l form ... does not satisfy the 
... conditions at t = 0, there being no derivative at that point ." Yet, worr ies about this fact prevented Amos ( 1969) from obtaining the solution 
in Eq. U4K though he had all the "pieces" of it. Meanwhile, Ting 1 1963) shows his solutions for flows of second-grade fluids have continuous 
derivatives at t = 0, however, all of these are for compatible initial and boundary conditions. 
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Oldroyd-B solution representation 1 ( IChristov and Jordan , 120091) . Define k :— A2I A\ and assume > 0. Using 
the Fourier sine transform in y, solving the resulting ordinary differential equation in t with the Laplace transform, 
one obtains 



erfc 



K< 1, 
K = 1, 



K> 1, 

(18) 



f * + Cf . sin {-§=) d^ + /| sin ( J=) df + £ + <£, sin ( J=) df 

0|erfc(^=), 

1-2JC" <Z/ + (^)sin(^=)d£ 

where < U ± (i;, t) correspond to fit;) ^ 0, respectively: 

_ exp[-^fM 1 ]{V7^cosh[(f//li) > /7©]+^sinh[(fMi) > //^]} exp[-g(£)t/Ai] sinh[(f/Ai) y7(?)] 



and 



exp[-*(0fMi] { VLf©icos[(/M!) VLf©l] + *<£) sinta/Aj) VI7X?)T]} *f rap[-*(£)tMi] sin[(t/Ai) VST^f] 



Oldroyd-B solution representation 2 (ITannei , Il962h . Define k :— A2IM and assume A12 > 0. Using the Laplace 
transform in t and the Bromwich integral inversion formula, one obtains 



u(y,t) = U H(t) 



where 



2 n Jo 



exp<- 



IvAa 



1 M(t])[ cos 6(j]) - sin 8(7])] 



x sin — 77 - 



/vA\ 



^M(rj)[cos0(r]) + sin6»(77)]l — I, (19) 



1 1 +T] 2 

1 + /c 2 ?7 2 ' 



6/(77) = 5 [tan 1 77 - tan '(707)] . 



Oldroyd-B solution representations (lMorrisonLll956l) . Define k :- A2/A1 and assume A\^ > and k < 1. Using the 
Laplace transform in t and a special splitting of the resulting quotient, one can invert the transform-domain solution 
using standard tables of inverses to obtain 



u(y, t) = U H(t) exp 



(1-*) t 



erfc 



2 Vvf V* 



i-K r tlM 

K JO 



2 VI 



V? 



2vr 



• exp 



(2-k) ( t 



erfc 



2 V V/ h<T V*. 



(20) 



where J p (-) is the Bessel function of the first kind of order p. 



Morrison! ( 1956 ) and Tannei ( 1962 ) both neglect to multiply their solutions by H(t). Though Morrison! ( 1956 ) was 



not studying the Oldroyd-B fluid specifically, he obtained Eq. (I17al i for the velocity in a viscoelastic rod whose stress 
res ponse is modeled b y a dashpot in series with an element consisting of another dashpot and a string in parallel. 



Vieru et al.l (120081) claim a "new" solution to the IB VP in Eq. (JT7J is obtained in the ir paper, though they use the 



Laplace transform and a quotient splittin g very similar to the one in dMorrison , 19561) . Additio nally, it is claime d 
that "the diagrams of the solutions" from (IFetecau and FetecaulEooll) 1 'are identical" to those in (IVieru et all 12008). 
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Though it remains unclear why a diagram (rather than an accurate plot of the solutio n) is relevant, the solution i n 



( I Vieru et all 120081) is suspect because it appears to agree with the incorrect solution in dFetecau and Fetecaul 120031) . 



The latter along with t he solution for the porous half-space and porous half-spa ce over a heat plate version s of this 
problem presented in dTan and Masuokal 2005blah were shown to be wrong by Christov and Jordan] ( 2009 ). Other 
studie s of Oldroyd-B dFetecaul 120021; iFetecau et all 120081: lKhanll2009HFetecau et al.[|2009h and Burgers Tkhan et al " 



201 Oh fluid flows are also erroneous because the mistake in applying the Fourier sine transform is made. 



3. Numerical solutions by finite-difference methods 

To provide an independent check on the transform solutions given in Sec. [2] we also solve the corresponding IB VPs 
numerically. The (uniform) spatial and temporal step sizes are defined as Ay := L/(M- 1) and Af := tf/{K- 1), where 
M > 2 and K > 2 are integers and now (y, t) e (0, L) x (0, t/\. Also, we let u" « u(jj, f) be the approximation to the 
exact solution on the grid, where yj := jAy (0 < j < M - 1) and f :- nAt (0 <n < K — 1). For appropriately chosen 
L » 1, the front does not reach the y — L boundary for any t e (0, ?/], and so this is the "numerical infinity." 

For the computations shown below, we use M — K — 5000 and L — 20 to obtain highly-accurate solutions. 
Matlab's built-in Gaussian elimination algorithm is used to invert the symmetric tridiagonal matrices resulting from 
the spatial discretizations. Additionally, the integral representations of the analytical solutions are evaluated using the 
high-precision numerical integration routine NIntegrate of the software package Mathematica (ver. 7.0.1). 

3.1. Second-grade fluid 



As shown in (IChristov and Christovl lin press!) , we may discretize Eq. ( 1 1 3at as follows: 

= vS n 8y- [± (uf 1 + u';)] + a6 y+ 6 t+ 6 y -u n j. 



(21) 



Here, 6 t + is the forward te mporal difference operator and 6 y+ and 6 Y - are, respectively, the forward and backward 
spatial difference operators ( IStrikwerdal |2004 §3.3). The boundary conditions from Eq. (1 1 3bb are implemented as 



[0, n = 0, 

\U , \<n<K-l; 



x M-\ 



0, 0<n<K-l. 



(22) 



The initial condition i s u° = 0(0< j < M — 1) owing to Eq. (1 1 3cb . It is a straightforward, though lengthy, calculation 
(see, e.g., Strikwerdal 2004 ) to show that this impli cit, two-level Crank-Nicolson-type discretization is uncondition- 
ally stable and has truncation error 0[( At) 2 + (Ay) 2 ] . lAmosI d 19691) gives some explicit schemes for Eq. d!3al ), however, 
the present one is superior in both its accuracy and stability. 

Since there exis t dimensionless variables t - t(v/a) and y = yj ^[a such that the problem and solution no longer 
depend on v and a ( Coleman et alin~965l : lBandelli et al. . 1995 ). the qualitative shape of the solution is invariant and 
varying the parameters is not enlightening in any way. Therefore, in Fig.[TJ we have made use of these dimensionless 
variables by showing u(y, T)/Uq rather than u(y, t). It is clear that the three analytical solutions to the IB VP in Eq. ( [TBI 
obtained by integral transforms presented in Sec. l2.1l agree identically with its numerical solution. 



3.2. Oldroyd-B fluid 

Constructing a reliable finite-difference scheme directly for Eq. (117at subject to Eq. ( 1 17bb turns out to be a difficult 
task. An easier approach is to introduce the fluid's acceleration w = du/dt with ro" ~ w(yj, f), then Eq. ( 1 1 Tab can be 
trivially rewritten as a system that we discretize by the following semi-implicit Crank-Nicolson-type procedure: 



X Y 8 t+ Xo) = -i(rof 1 + »») + vVM- + " 2 + vAaVV [\ {^f + »")] • (23) 
Noting that w(0, t) = Uo6(t), the boundary conditions from Eq. d!7bb are implemented as 

n = 0, 

«Sf-i = »2r-i = 0, 0<n<K-l. (24) 



,n-l/2 



0, n = 0, 



\U , \<n<K-\; 



4tt, 1 < n < n\ 

n*At ' ' 

n* <n<K- 1; 
6 



u/Uq 




Figure 1: Solutions to Stokes' first problem for the second-grade fluid for t = 1 (bottom curve) and t = 10 (top curve). Legend: Fourier sine 
transform solution from Eq. j!4t (•), Puri's Laplace transform solution from Eq. )15t (■), Bandelli et al.'s Laplace transform solution from 
Eq. )16t (♦) and the numerical solution using the scheme from Eq. )2U (-). 



The initial conditions are uj 1 ^ 2 = ro'! = (0 < j < M — 1) owing to Eq. d 17cb . The scheme is not sensitive to the 
parameter n* ', so we take n* = 10 in our calculations. 

Unfortunately, due to the implementation of the (5-function boundary condition, it is no longer straightforward to 
show stability. Nevertheless, numerical experiments show the scheme is stable for At = O(Ay). It is easy to establish 
the truncation error is 0[(At) z + (Ay) 2 l. lTownsendld 19731) gives a 0[(At) 2 + Ay] fully-implicit scheme for the system 
of Eqs. ( [Tol l and (fTTT i. however, the present scheme is si mpler and more accurate. Oth er modern numerical approaches 
to one-dimensional viscoelastic flows can be found in ( Amoreira and Oliveira . 2010b an d the references therein. 

By using the dimensionless variables t = t/A\ and y — yl yjvA\ ( Tannerl 1962 ). it can be shown that quali- 
tative differences in the shap e of the solution result only from having k < 1, k = 1 or k > 1. Experiments 
(IToms and Strawbridgel 119531) suggest that k{= At/Ax) is in the range 0.05 to 0.4. Therefore, in Fig. [2] we take 
k = 0.2 and use the dimensionless variables above to show that the three analytical solutions to the IB VP in Eq. ( fTTl i 
obtained by integral transforms agree identically with its numerical solution. 




Figure 2: Solutions to Stokes' first problem for the Oldroyd-B fluid with k = 0.2 for t = 1 (bottom curve) and t = 7 (top curve). Legend: Fourier 
sine transform solution from Eq. 1181 (•), Tanner's Laplace transform solution from Eq. 1191 (■), Morrison's Laplace transform solution from 
Eq. (20} (♦) and the numerical solution using the scheme from Eq. (23) (-). 



4. Conclusion 



In the last decade, a disturbing trend has emerged in non-Newtonian fluid mechanics. Many authors have been re- 
deriving known results, specifically those from the 1950s through 1980s on the simple flows of certain non-Newtonian 
fluids. Unfortunately, many errors have been made in these "new" derivations. When the problems considered are 
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actually novel, it is usually a minute change (e.g., in a boundary condition) that distinguishes them from the classical 
works. Even when free of mathematical errors, these make little (if any) contribution to the mechanics of fluids. 

Here, we presented three known correct solutions (by the Fourier sine transform, by the Laplace transform with 
the Bromwich integral inversion formula and by Laplace transform with quotient splitting and tables of inverses) to 
Stokes' first problem for the second-grade and Oldroyd-B fluids. Additionally, we presented a representative list of 
the papers in which the so-called "new solutions" are wrong. A complete list would be far too long to attempt here. 
The numerical schemes constructed in Sec.[3]provide a wholly independent check on the classical Laplace transform 
solutions and the corrected Fourier sine transform solutions. 

Given the astonishing amount of misinformation on this topic in the literature, we offer some advice to future 
researchers based on the present work: 



The Laplace transform always works on a well-posed linear IB VP. The solution of lPuril (119841) remains the first 
correct exact solution to Stokes' first problem for the second-grade fluid, despite the denigrating remarks many 
authors, whose work is erroneous, make about this solution. 

-it 

The assumption that the fluid is initially at rest dictates that u(y, 0) = (identically), whence 'jp?(y, 0) = for 
any k. Using the latter equ ality for k > 1 does not constitute an "addit ional" or "unphysical" assumption as 
incorrectly claimed in (e.g.. flan and MasuokalEoolbT: IVieru et all 120081) and many derivative works thereof. 



• For partial differential equations with mixed derivatives and incompatible initial and boundary data, the solution 
does not have to satisfy the initial condition in backward time, i.e., as t — > + . 

• One has to be very careful in applying the Fourier sine (or cosine) transform to partial differential equations 
with mixed derivatives and incompatible initial and boundary data because distributional derivatives of the 
Heaviside function will inevitably have to be taken for Stok es' problems. Note that the same error exposed in 
dChristov and Jordan , 120091: IChristov and ChristovL lin press!) is committed when solving the unsteady version 
of Stokes' second problem for non-Newtonian fluids. However, when U(t) ~ r as t — > (k > 1), e.g., 
U(t) = sin(o;fX the singularity is ameliorated thereby allowing an erroneous derivation to produce a correct 
solution. The same is true for the "constantly accelerating plate" problem in which U(t) = At. 

• The decay boundary condition that u — » "sufficiently fast" as y — > oo is typically enough to guarantee a 
solution by the Laplace transform. In the case of the Fourier sine transform, keeping in mind we seek a classical 
solution of the partial differential equation, the implicit assumptions are made that u is twice continuously 
differentiable in y, and that u, du/dy and d 2 u/dy 2 are all integrable for y € (0, oo). (Note that this presupposes 
nothing about the regularity of u in f, which we saw above is quite lo w.) In fact, it is known from the theory of 
integration (see Theorem 33.7 and Proposition 23.8 in lPriestlevl 1 1 9971) that these (implicit) assumptions not only 
guarantee that the Fourier sine transform can be applied but also that u and du/dy — » as y — > oo. Confusion 
about this has led some authors to impose the additional (unnecessary) condition that du/dy — » as y — > oo. 

• Taking the limit as the non-Newtonian parameter(s) go to zero should reduce any solution to the known New- 
tonian one. However, this a necessary but not sufficient condition. Consequently, this exercise provides any 
insight only if one fails to recover the Newtonian solution, meaning the non-Newtonian one is wrong. For 
the problems considered in the present work, it happens that both the correct and erroneous non-Newtonian 
solutions reduce to the correct Newtonian one, which is unaffected by the error in computing fj(0, t) as its gov- 
erning equation does not have a mixed third derivative, hence one never has to compute ^-(0, t) when applying 
the Fourier sine transform. 

• The signs, thermodynamic restrictions on, and orders of magnitude of the non-Newtonian parameters cannot 
be ignored. Introducing dimensionless variables is an easy remedy to the lack of precise measurements of, e.g., 
the second-grade fluid's parameter a\ . Then, plots can be made with ease and without loss of generality. 

Another important class of non-Newtonian fluids are those of the Maxwell type, which exhibit stress relaxation 
but have no strain memory (ret ardation), i.e. , A2 = in Eq. A thorough overview of the various (correct) solu- 
tion representations is given by IJordan et al.l (120041) . And, a positivity -preserving numerical scheme, which can be 
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us ed as an independent check on any suspicious "new" solutions that may appear in the literature, was constructed 
bv lMickens and Jordanl ( 120041) . Another approach to ascertaining the correctness of vari ous solutions to Stokes' first 
probl em is to consider the asymptotic searin gs of the velocity and shear stress with time ( Muzvchka and Yovanovichl 
2010h that they predict. Finally, we note that lPreziosi and Joseph! dl987h also provide correct solu tions to Stokes' first 



proble m for viscoelastic fluids with a variety of memory kernels using the Laplace transform, while lPhan-Thien and C hew 
d 1 988b present numerical solutions to Stokes' first problem for a class of viscoelastic fluids that reduce to the Oldroyd- 
B and Maxwell models in certain distinguished limits. 
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